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CHAPTER  1 
INTRODUCTION 


bMI 


S' 


The  principal  objectives  of  this  program  were:  (1)  to  develop  a 
technique  for  the  measurement  of  isochoric  press ure~energy  relations  for 
solid  and  porous  materials,  and  (2)  to  exercise  this  technique  on  porous 
materials,  including  a metal  with  two  initial  densities  and  two  cermets. 

The  former  objective  is  the  subject  of  Part  I of  this  report;  the  porous 
material  results  are  presented  and  discussed  in  Part  II. 

For  instantaneous  energy  deposition  in  a thermoelastic  solid  the 
pressure-energy  relationship  is  linear  and  is  defined  by  the  relation 

P = YpE  (1) 

where  ? is  the  pressure  increase,  E is  the  internal  energy  change,  p is 
the  density,  and  Y is  the  Gruneisen  parameter.  Actual  energy  depositions 
I require  a finite  time  which  results  in  some  stress  relief  and  expansion 

t near  boundaries;  only  internal  material  may  have  been  heated  at  constant 

I voliB&e.  For  composite  and  porous  materials  the  pressure-energy  relation 

is  non-linear,  and  the  above  expression  defines  an  "effective”  Gruneisen 
! paramcirer  wtiich  is  a fvmctlon  of  energy  and  Initial  volume.  As  the  intcnml 

energy,  E,  is  Increased  at  constant  volume  in  a porous  taatorial,  the 
pressure,  P,  initially  increases  somewhat  proportionally  to  E;  at  higher 
energy  the  matrix  material  yields  and  expands  into  the  voids  and  there 
is  no  further  pressure  Increase  (until  voids  arc  filled).  At  oven  higiter 
energy,  the  vapor  state  of  the  taatrlx  is  reached  (or  the  voids  become 
filled)  and  the  pressure  increases  rapidly  with  energy.  This  complex  P-6 
behavior  is  illustrated  in  Figure  1.  Also  shovti  in  Figure  1 is  the  linear 
behavior  for  a selected  thermoelastin  solid  tdrosc  slope  (Gruneisen)  is  less 
; than  the  initial  slope  for  the  porous  material.  At  some  point  (i.e.,  where 

the  curves  cross)  the  solid  and  porous  materials  develop  the  same  pressure 
I for  tlte  same  change  in  internal  energy.  Thus  a judicial  cltoice  of  solid 

< material  and  internal  energy  permits  matching  the  porous  material  stress. 
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Energy 


Figure  1«  Preiisure~Caergy  Behavior  for  Solid  and  Porous  Materials 


Also  of  importance  is  the  adiabat  centered  at  the  initial  pressure- 
energy  state  achieved  by  instantaneous  energy  deposition  at  constant  volume, 
as  indicated  in  Figure  2.  In  general,  because  Gruneisen  measurement  experi- 
ments involve  stress  wave  measurements,  the  adiabats  are  inextricably 
related  to  the  determination  of  the  initial  pressure-energy  or  Gruneisen 
state  as  discussed  below. 

Most  p.:essure-energy  or  Gruneisen  parameter  determinations  by 
electron  beam  (e-beam)  deposition  have  utilized  relatively  thick  samples; 
i.e.,  thick  with  respect  to  the  electron  range.  An  example  of  this  technique 
is  described  by  reference  1.  The  energy  gradient  in  the  sample  and  the 
finite  deposition  time  result  in  a poorly  defined  "initial"  pressure-energy 
state,  and  in  most  cases  only  a small  region  of  the  sample  is  heated  at 
constant  volume.  A complex  material  equation  of  state,  a porous  material 
model,  and  a wave  propagation  computer  code  are  generally  required  in  order 
to  compute  the  "initial"  pressure-energy  distribution  from  tJie  measured 
response.  The  stress  attenuation  which  occurs  during  wave  propagation 
between  the  heated  portion  of  the  sample  through  tlui  cold  sample  to  the 
stress  or  particle  velocity  gauge  at  the  rear  of  the  sample  limits  the 
resolution  of  the  initial  sample  pressure  and  hence  the  determination  of 
the  Gruneisen  parameter. 

In  more  recent  measurements  this  situation  has  been  considerably 
improved  through  the  use  of  uniformly  heated  sampler  utilising  relatively 
high  mean  electron  energy  o-beams  and  fused  silica  buffers  (see  for  example, 
reference  2).  Tlie  extremely  low  fused  silica  Gruneisen  parameter  results  in 
relatively  low  stress  generation  in  that  material  which  transmits  tlte  stress 
wave,  tite  fimed  silica  velocity-time  profile  is  measured  wtrh  an  inter- 
ferometer, or  the  stress  wave  is  measured  with  a quartz  gau^e.  The  knowl- 
edge of  the  silica  adiabats  allows  a point  on  the  sai^le  expansion  adiabat 
to  be  directly  determined  without  resorting  to  complex  computer  computations; 
but  the  adiabatic  path  joining  this  point  to  the  initial  state  is  generally 
not  known  and  must  he  assumed  in  order  to  find  the  initial  (Gruneisen) 
pressure.  The  determination  of  the  internal  energy  change  (dose)  is  also 
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Figure  2.  Pressure'-Volume  Adiabat  from  Initial  Constant  Volume  State 


generally  not  accurately  known  because  it  depends  on  the  electron-beam 
fluence  and  depth-dose  profile,  which  in  turn  depend  on  calorimetric 
measurements  and  machine  repeatability. 

In  the  technique  utilized  in  this  work,  a stress  generating 

reference  material  replaces  the  fused  silica,  and  fairly  uniform  and  rapid 

deposition  is  desired  through  both  sample  and  reference  material,  as 

indicated  in  Figure  3.  This  results  in  a significant  improvement  in  the 

ability  to  determine  the  initial  pressure-energy  state.  The  measured 

particle  velocity-time  history  at  the  rear  surface  of  the  reference 

material  is  characterized  by  two  "steps"  of  total  amplitudes,  U and  U , 

t s 

as  indicated  in  Figure  3.  Simple  analysis,  presented  in  Section  5, 
allows  a point  on  the  sample  adiabat  close  to  the  Initial  (Gruneisen) 
pressure-energy  state  to  be  determined  from  U and  U . Tlie  first,  or 
reference  step,  amplitude  determines  the  internal  energy  change*,  and 
the  second,  or  sample  step,  amplitude  determines  the  adiabat  point. 

Several  points  on  an  adiabat  (both  compression  and  expansion)  can  be 
found  by  utilizing  various  reference  materials  while  maintaining  the  same 
e-beam  dose. 

The  sample  pressure  is  found  by  a differential  technique  which 
compares  the  reference  and  sample  step  amplitudes.  In  the  optimal 
Gruneisen  experiment  the  sample  steps  are  of  equal  amplitude  (i.e., 

U • U ),  implying  that  the  sample  Gruneisen  stress  equals  the  known 

r 3 

reference  material  Gruneisen  stress.  This  is  optimal  because  when  the 
reference  material  Gruneisen  stress  equals  the  sample  Gruneisen  stress 
no  extrapolation  along  the  sample  adiabat  is  required.  It  is  clearly 
desirable  to  design  the  experiments  so  as  to  minimize  the  stress  dif- 
ferential between  sample  and  reference  materials.  This  new  technique 
will  be  referred  to  as  the  "stress-matching"  technique. 


* It  was  shown  in  1967  (reference  3)  that  calorimetric  e-beam  measurements 
could  be  made  via  stress  measurements  in  thermoelastic  absorbers. 
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CHAPTER  2 
EXPERIMENT  DESIGN 

In  order  to  optimally  exercise  the  "stress-matching"  technique 
Introduced  above  one  must  select  appropriate  reference  materials  and  the 
e-bean  machine.  These  selections  are  coupled  to  the  thermomechanical 
properti'^s  of  eich  particular  sample  material;  i.e.,  the  same  reference 
materials  and  electron  beam  will  not  provide  the  optimum  experiment  for 
all  porous  sample  materials  and  dose  ranges.  However,  some  general  guide- 
li’'es  and  procedures  for  the  design  of  optimal  experiments  can  be  given. 

The  principal  requirement  for  reference  materials  are:  (1)  that 
the  material  respond  in  a known  and  reproducible  manner,  and  (2)  that  the 
material's  Gruneisen  pressure  gs.neration  (product  of  Gruneisen  parameter 
and  density)  for  a gj.«/en  internal  energy  increase  (dose)  is  close  to  that 
of  the  sample  to  provide  a near  stress-matching  condition.  This  implies 
that  a number  of  refe^'ence  materials  are  required  to  characterize  a single 
porous  sample  material  over  an  energy  range.  r.»e  general  validation  of  the 
former  reauiromcnt  for  u number  of  reference  materials  is  the  subject  of 
Section  6.  Candidate  reference  materials  vsre  selected  to  provide  a range 
of  prcssure^energy  reeponses  in  order  to  satisfy  the  latter  requirement. 

Figure  4 indicates  the  press uve-energy  charactarlstics  of  the 
selected  reference  materials  and  the  range  of  pre8su;.&*energy  response 
anticipated  for  the  porous  sample  materials  (shaded  region).  Tlie  sample 
pressure  will  bo  matched  by  the  various  reference  nuiterials  at  the  dose 
levels  corresponding  to  the  abscissas  at  the  intersections  of  the  reference 
and  sample  response.  The  reference  materials  selected  were:  (1)  aluminum 
(6061-T6),  (2)  quartz  (X-cut),  (3)  germanium  till),  (4)  silicon  (poly- 
crystaline),  and  (5)  silica  (fused). 

Tt-^e  e-beom  machine  choice  is  based  upon  meeting  the  requirements 
for  spatial ly-vide  uniform  pressure  regions  in  the  samp'^e  and  reference 
materials,  the  response  of  which  produces  a well-oefined  particle  veloc'.ty 
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Energy,  E (cal/g) 

Figure  4.  NcHQinal  Pressure-Energy  Characteristics  for  Selected 
Reference  Materials  and  the  Anticipated  Sample 
Characteristics 
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history  (stress  profile).  The  electron  beam  depth-dose  profile  and  the 
deposition  time  both  control  these  requirements. 

A simple  computer  routine  was  xised  to  predict  the  energy,  stress, 
and  particle  velocity  profiles  from  several  electron-beam  machines  for 
the  various  reference  materials.  The  inputs  required  are; 

1.  depth-dose  profile 

2.  deposition  time 

3.  sample  and  reference  material  density,  wave  speed, 

Grunelsen  parameter,  and  thickness 

Figures  5 and  6 are  examples  of  the  predicted  stress  profiles  which  exist 
in  the  sample  and  reference  material  just  after  deposition.  Figure  5 
indicates  the  dependence  on  mean  electron  energy  (deposition  profile) . The 
8 MeV  Hermes  beam  provides  the  required  regions  of  uniform  pressure  while 
the  lower  4 MeV  FX-75  beam  results  in  a triangular  stress  wave,  the  amplitude 
of  which  depends  on  the  time  dependent  (rather  than  time  integrated)  depth- 
dose  profile  and  on  stress-wave  attenuation.  Figure  6 contrasts  the  results 
for  the  Hermes  II  and  the  Aurora  beams,  both  nominally  8 MeV  beams  with  70 
and  140  ns  (FWHM)  deposition  times,  respectively.  These  examples  used  the 
following  representative  sample  properties:  density  - 14  g/cm  , wave  speed  - 
0.5  cm/ps,  Grunelsen  parameter  - 0.3.  X-cut  quartz  was  the  reference  material. 

Based  upon  these  types  of  considerations,  the  Hermes  IX  accelerator 
was  selected  for  tlie  Grunelsen  experiments. 


0 0.1  0.2  0.3  O.A 

Time  (ps) 

Figure  S.  Predicted  Response  of  a Typical  Target  Configuration 
Indicating  Dependence  on  Mean  Electron  Energy 
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Figure  6 


Predicted  Response  of  a Typical  Target  Configuration 
Indicating  Dependence  on  Deposition  Time 


CHAPTER  3 

EXPERIMENT  CONFIGURATION 

i 

This  section  describes  the  physical  layout  of  the  experiment 
apparatus,  and  presents  examples  of  the  Interferometer  data  which  were 
recorded  for  the  reference  materials. 

The  target  configuration  used  at  the  Sandia  Laboratories*  Hermes  II 
accelerator  is  shown  in  Figure  7.  Either  a sample  holding  cup  or  a five- 
element  spatially-resolving  graphite  calorimeter  was  positioned  in  the 
center  of  the  target  assembly.  The  outer  annular  area  contained  twelve 
calorimeter  elements  positioned  on  three  radii  as  shewn.  The  calorimeter 
thermocouple  output  voltage  was  recorded  and  reduced  by  the  Sandia  computer 
linked  data  acquisition  system  (DADS),  reference  4.  Representative  fluence- 
radius  data  are  presented  in  Section  4.  The  velocity  of  the  free  rear 
surface  of  the  reference  material  was  measured  with  a cemabination  displace- 
ment (reference  5)  and  velocity  laser  interferometer  (reference  6),  The 
velocities  in  the  range  below  0.1  mm/ys  were  best  resolved  by  the  displacement 
interferometer  and  higher  velocities  with  the  velocity  interferometer. 

Greater  than  normal  velocity  resolution  and  the  ability  to  resolve  otherwise 
subtle  velocity  reversals  were  obtained  with  a quadrature  system  which  pro- 
vides two  phase  shifted  Interferometric  signals.  These  signals  were 
transported  about  50  feet  by  a single  laser  beam  from  the  Interferometer  to 
a shielded  room  where  they  were  optically  separated,  detected  with  pulsed 

j photomultipliers,  and  recorded  oscillographlcally.  A second  laser  beam 

transported  the  displacement  interferometer  signal  to  the  shielded  room. 

! This  overall  setup  in  the  Hermes  II  facility  Is  shown  in  Figure  8.  The 

oscillograph  and  the  resultant  velocity-time  history  shown  in  Figure  9 are 
representative  of  the  velocity  Interferometer  data.  The  points  in  the  Figure 
correspond  to  the  location  of  the  intensity  maxima  nnd  minima  for  one  quad- 

I 

I rature  signal.  Tlte  velocity  increment  between  adjacent  intensity  extrema 

I is  a constant  given  by  (reference  6) 

I - If  (2) 
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Figure  7.  Target  and  Calorimeter  Configuration 
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Figure  8.  Set  Up  at  ttie  Hermes  II  Electron  Beam  Accelerator 
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uhere  X and  c are  the  wavelength  and  speed  of  light  and  L Is  the  length 
of  an  optical  delay  path.  The  velocity- time  relation  between  maxima  and 
minima  was  obtained  from  the  analytic  relation  between  velocity  and 
intensity  (reference  7). 

U(t)  { cos"^[2I(t)  - 1]}  (3) 

The  velocity  data  shown  in  Figures  10  and  11  were  obtained  by 
differentiating  the  displacement  interferometer  data.  The  points  in  the 
figure  represent  the  average  velocity  over  the  time  interval  between  N 
adjacent  intensity  extrema»  which  corresponds  to  a displacement  of 

Ax  » NX/4  (4) 


I 

1 


[. 


The  velocity  is,  dierefore, 


^ 

At  " 4tAt 


(5) 


In  the  three  above  representative  examples  tite  velocity  begins 
to  Increase  as  the  beam  begins  to  deposit  energy,  and  the  reference 
plateau  is  reached  about  100  ns  later  when  deposition  is  nearly  complete. 
The  second  step  is  the  result  of  the  arrival  of  the  stress  wave  from  the 
sample. 

The  characteristic  "reference"  and  "sample"  material  particle 
velocities  (U  and  U„)  are  Indicated  on  the  figures,  these  velocities  are 
selected  and  used  in  the  analysis  procedure  as  described  in  Section  S to 
determine  the  inteneediate  adiabatic  expansion  stress  state,  which  is 
close  to  the  Gruneisea  stress,  and  thence  tite  Grunelsen  stress  state  of 
interest. 
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Figure  10.  Representative  Dlsplaceaeut  Interferotteter  Oats  for 
Qtz/Si  Ck>n£iguratioo  - Shot 
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CHAPTKR  4 

ELECTRON  BEAM  CHARACTERIZATION 

The  He:  mes  II  pulsed  electron-be/'in  machine  was  used  as  the  electron 
source  for  all  the  experiments  described  in  this  report.  The  machine  and 
facilities  are  described  in  detail  in  reference  8 and  the  general  performance 
data  are  «ivan  iu  Table  1.  Diode  current  and  voltage  data  were  taken  for 
each  of  about  100  shots  conducted  under  this  program. 

The  depth-dose  profile  was  determined  from  the  thermoelastic  response 
of  a 0.50-inch  thick  X-cv\t  quartz  crystal.  The  rear  free  surface  particle 
velocity  time  history  was  measured  with  a laser  velocity  interferometer. 
Transport  calculations  were  made  with  the  ELTRAN  (reference  9)  Monte  Carlo 
code  using  the  electron  energy  spectrum  determined  from  the  diode  current 
and  voltage  diagnostic  data  for  the  shot  and  various  Incident  angles.  The 
measured  and  computed  depth-dose  relations  are  normalized  to  the  absorbed 
energy  and  compared  in  Figure  12.  The  best  fit  to  the  experimental  data 
was  obtained  for  an  assximed  incidence  angle  of  30  degrees.  This  angle  was 
used  for  all  other  calculations  presented  in  this  report. 

Calculations  were  made  with  several  representative  spectra  determined 
from  diode  data  representing  the  range  of  mean  electron  energies  obtained 
throughout  the  tests,  and  the  peak  dose  (energy  coupling  coefficient)  was 
observed  to  vary  by  i 5 percent. 

The  spectrum  shown  in  Figure  13  for  shot  No.  139  was  chosen  as 
producing  the  most  representative  depth-dose  profile  and  all  calculations 
utilize  this  spectrum.  Also,  all  other  calculations  are  normalized  to 
incident  fluence  and  all  finance  data  arc  corrected  to  incident  fluence  as 
described  below.  Figure  14  presents  the  ELTRAN  results  for  this  spectnua 
incident  dt  30  degrees  to  darbon  (atomic  number,  Z ■ 6),  germanium  (Z  «•  32) 
and  tantalum  (Z  73).  Calculations  for  actual  bi- laminate  target  con- 
flguratjans  are  presented  in  Sections  6 for  the  reference  material 
experiments  and  in  Part  XI  for  the  porous  sample  configurations.  The 
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TABLE  1 


HERMES  II  ACCELERATOR  GENERAL  CHARACTERISTICS 


Marx  Generator  Voltage 

70.0  kV 

Peak  Diode  Voltage 

9.4  MV 

Peak  Diode  Cxirrent 

11.5  U 

Mean  Electron  Energy 

7.8  MeV 

Pulse  Width  (FWTO!,  power) 

70  ns 

Total  Beam  Energy 

76  kj 

Normalized  Dose,  E/(J) 


2 

Depth  (g/an  ) 


Figure  12.  Measured  and  Computed  Depth-Dose  Relations  for  Quartz 


2 


3 


Depth  (g/m  ) 


Figure  14.  Depth-Dose  Calculations  for  Three  Materials 


deposition  profile  was  flat  to  within  5 percent  across  the  sample,  with 

2 

the  aid  of  a thin  0.176  g/cm  carbon  filter. 

The  deposited  energy-versus-time  characteristics  as  determine  i 
from  diode  data  are  shown  in  Figure  15.  The  70  ns  linear  approximation 
shown  is  the  full  width  at  half  maximum  of  the  power  curve,  and  is 
suggested  for  use  in  time-dependent  material  response  calculations. 

Fluence  data  were  taken  with  a spatially  resolving  graphite 

calorimeter  (shown  previously  in  Figure  7) . The  center  five  elements  were 

interchangeable  with  the  sample  holding  cup  so  that  peripheral  fluence 

measurements  could  be  made  during  sample  shots.  All  the  absorbed  fluence 

data  reported  here  are  converted  to  incident  fluence  based  on  electron 

deposition  calculations  for  graphite  which  show  that  a total  of  12.4  percent 

of  the  incident  energy  is  lost  by  absorption  in  the  0.1  cm  thick  filter,  by 

Bremsstrahlung  and  by  backscattering.  Three  representative  sets  of  fluence- 

radius  data  are  given  in  Figure  16  for  shots  with  peak  fluences  of  300,  150, 
2 

and  50  cal/cm  . 

The  various  fluence  levels  were  obtained  by  adjusting  the  axial 

distance  from  the  anode  to  the  target/calorimeter  assembly.  The  highest 

2 

fluence,  about  300  cal/cm  , was  obtained  at  the  "first  pinch",  which  occurs 
at  about  13  cm  from  the  anode  plane.  This  is  less  than  half  of  the  peak 
fluence  anticipated  based  on  data  presented  in  reference  10,  which  were 
taken  with  a larger  diameter  drift  tube  (71  cm  versus  our  18  cm  tube), 
and  this  may  have  resulted  in  a more  intense  pinch.  Thus,  the  peak  dose 
in  our  experiments  was  about  100  cal/g. 
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CHAPTER  5 

ANALYSIS  PROCEDURES 


Tills  Section  describes  the  procedures  used  to  compute  the  initial 
sample  pressure-energy  state  and  other  related  parameters  of  interest. 

The  rear  free  surface  velocites  and  U„  are  defined  by  the 

amplitude  of  the  first  and  second  velocity  steps  respectively.  Since 

the  deposition  is  not  instantaneous  and  uniform  through  the  thickness  of 

the  sample  and  reference  materials,  the  step  amlitudes  are  determined  as 

indicated  in  Figures  9,  10,  and  11.  The  maximum  velocity  occurring  near 

the  center  of  the  "sample  step"  is  taken  as  U . The  velocity  found  by 

extrapolating  the  reference  step  to  the  sample-reference  interface  is 

taken  as  U„. 

r 

Figure  17  indicates  the  adiabatic  paths  which  relate  the  measured 
step  amplitudes  U_  and  U„  to  the  initial  (Gruneisen)  pressures  P and  P„. 

The  sample  adiabat  is  represented  by  a solid  line  and  the  reference  material 
adlabats  by  dashed  lines.  All  reference  material  adiabats  are  assumed  to 
be  linear  for  purposes  of  this  discussion.  The  sample  adiabat  is  assumed 
to  be  linear  between  P and  P„.  Note  the  existence  of  intermediate  states 
which  are  defined  by  the  Hugoniot  continuity  and  conservation  requirements. 
The  intermediate  state  on  the  sample  adiabat  (U„,  P ) can  be  found  from  the 
measured  velocities  U and  U.  and  the  known  reference  material  adiabats. 

fc  8 

For  the  case  of  linear  reference  material  adiabats  with  constant  slope 
the  relations  are  simply 


Pa  “ 1/2  Zy{Uy  + Ug)  (6) 

and 

“a  - 1/2  ("a  - "j)  O) 

Tltese  relations  define  a point  on  the  sample  adiabat  centered  at  the 
initial  heated  state.  Note  that  various  points  on  the  sample  adiabat  from 
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the  same  initial  P_,E  state  may  be  determined  by  repeating  the  experiment 

with  different  reference  materials  which  generate  different  initial 

pressures  Pj.*8»  and  have  different  adiabat  slopes.  If  the  initial  reference 

material  pressure  is  greater  than  the  initial  sample  material  pressure 

(P_  > P_),  a point  on  the  sample's  compression  curve  is  determined.  In 

this  case  the  reference  material  is  placed  in  tension,  vdilch  may  result  in 

fracture  (spall)  of  the  reference  material,  precluding  the  measurement  of  U at 

8 

the  reference  rear  surface.  Figure  11  is  an  example  of  such  a case  which 
did  not  result  in  spallation.  If  the  Initial  reference  material  pressure 
is  less  than  the  initial  sample  material  pressure  (P  < P ),  a point  on 

T 8 

the  expansion  adiabat  is  determined.  By  using  different  reference  materials, 
a complete  adiabat  can  be  mapped  from  the  initial  P-,E^  state. 

8 S 

If  the  sample  adiabat  shape  between  P and  P is  known,  the  initial 

3 S 

sample  pressure  P„  can  also  be  obtained  from  U and  U . If  the  portion 
of  the  sample  adiabat  connecting  the  initial  state  with  the  intermediate 
state  is  assumed  to  be  linear  with  slope  in  the  P - U plane,  the  initial 
sample  (Gruneisen)  pressure  is: 

Pg  - Py  + i/2CZj.  + ZgXUg  - U^)  (8) 

where 


P - Z U 
r r r 


(9) 


For  linear  elastic  san^le  and  reference  materials,  the  adiabat 
slopes  are  simply 


Z - p D 
8 *^8  8 


(10) 


and 


Z p D 
r '^r  r 


(11) 


where  p and  0 ere  the  Initial  density  and  the  dilatatlonal  wave  velocity 


in  the  direction  of  wave  propagation.  In  order  to  simplify  data  analysis, 
it  is  highly  desirable  to  utilize  linear  elastic  reference  materials. 

For  poroxis  or  composite  sample  materials  the  adlabat  slopes  are  not 
constant  and  the  local  wave  speed  can  vary  more  than  a factor  of  two  between 
the  initial  P ,E  state  and  the  expanded  zero  pressure  state.  In  general, 
the  slope  along  the  adlabat  is  given  by 

Zs(?)  - § > P3IXP)  (12) 


\diere  D(P)  is  the  local  Lagrangian  wave  speed.  E-beam  wave  transit  time 

data  were  used  to  obtain  an  effective  adlabat  slope  Z for  the  segment 

s 

between  P and  P as  follows, 
a s 

After  deposition,  relief  waves  propagate  toward  each  other  from  the 

sample  material  surfaces  meeting  in  the  center.  The  velocity  of  the  leading 

edge  of  these  waves  is  D(P  ) . From  the  center  of  the  sample  a wave  then 

6 

propagates  to  the  elastic  reference  material  at  velocity  D(P,).  Tlie  average 

o 

wave  speed  found  by  measurement  of  the  sample-step  pulse-width  is  therefore 
(neglecting  the  very  small  particle  velocities): 


D(P  ) + D(P  ) 
X _ ^ 8 ^ a^ 

T 


(13) 


where  x is  the  sample  thickness.  Therefore,  D confuted  as  x/t  may  be 

s 

used  to  obtain  an  effective  adlabat  slope  for  use  in  the  linear  analysis 
scheme.  Equation  (8).  11:g  experimental  relief  wave  speeds  agreed  with 

the  literature  values  (see  Chapter  6,  Table  3). 

As  the  pressure  gotiorated  in  the  reference  taaterial  P^.  approaches 

the  sample  pressure  P the  portion  of  the  sample  adiabat  joining  the 

s 

desired  P„,E„  state  with  the  point  on  the  sample  adiabat  becomes  shorter 

S 8 

and,  therefore,  titc  uncertainty  in  the  measurement  of  P decreases 

9 

(dependence  on  Z ■ pD  is  relaxed.)  It  is  therefore  desirable  to  use  a 
8 8 

reference  material  which  closely  matciws  the  pressure  generated  by  the  sample. 
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The  change  in  Internal  energy  for  a particular  test  can  be  directly 
determined  from  the  initial  reference  material  response  since 

P 

(14) 

and  is  known  by  Equation  (9) . If  the  sample  and  reference  materials 
achieve  slightly  different  energy  states  because  of  non-uniform  deposition, 
the  sample  energy  state  is  determined  from  the  coupling  coefficient  ratio 

r 

The  coupling  coefficelnts  are  defined  as  follows: 

Cs  ■ V*1 

and 


Cy  - (17) 

where  is  the  incident  fluence.  'Die  appropriate  coupling  coefficients 
are  found  from  ELTRAM  calculations  (Section  4)  and  are  selected  to 
correspond  to  the  measured  particle  velocities  and  U„,  as  described 
in  Section  3i  i>e.,  C is  the  peak  value  and  C is  the  value  at  the  sample- 

a ft  . 

reference  interface.  (Ttte  values  for  the  reference  material  test 
configurations  are  given  in  Section  6,  table  2.) 
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CHAPTER  6 

RESULTS  FOR  REFERENCE  MATERIALS 


In  order  to  exercise  the  stress-matching  technique  on  porous  and 
composite  materials  of  interest,  the  thei.aomechanical  response  of  a group 
of  appropriate  reference  materials  was  verified.  Several  candidate 
reference  materials  were  chosen  based  on  the  anticipated  pressure-energy 
response  of  the  porous  sample  materials,  the  response  of  which  is  dis- 
cussed in  Part  II.  These  reference  materials  were  evaluated  (calibrated) 
by  performing  therracmechanical  response  measurements  in  the  previously 
described  test  configuration  in  which  two  candidate  reference  materials 
were  exposed  in  each  test.  One  material  in  the  "sample"  position  (nearest 
to  the  pulsed  electron  source)  and  the  other  in  the  "reference"  position. 
The  rear  free  surface  particle  velocity  of  the  reference  material  was 
measured  to  provide  the  step  particle  velocity  amplitudes  U and  U as 
described  previously. 

The  energy  coupling  coefficients  were  determined  from  ELTRAN  code 
Monte  Carlo  deposition  calculations  based  on  beam  characteristics 
established  in  Section  4.  Deposition  calculations  were  made  for  each  test 
configuration.  Representative  deposition  profiles  are  given  and  compared 
with  the  nominal  carbon  deposition  profile  in  Figures  18  to  20.  Table  2 
summarizes  the  deposition  profile  data  in  terms  of  the  energy  coupling 
coefficients  and  for  each  test  configuration. 

Based  on  the  analysis  scheme  described  in  Section  5,  the  various 
pressure,  energy,  and  particle  velocity  states  were  computed  from  the 
coupling  coefficients  and  measured  free  surface  particle  velocities.  Die 
reference  material  parameters  which  were  used  in  the  data  analysis  are 
given  in  Table  3. 

The  measured  and  computed  data  for  the  reference  materials  are 
summarized  in  Table  4,  and  the  pressure-energy  data  are  plotted  in 
Figure  21.  The  Grunelsen  parameters  were  obtained  from  the  slope  of 
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Normalized  Dose,  E/<|> 


Figure  20,  Deposition  Profile  for  Ge/Qtz  Configuration 
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TABLE  2 


ENERGY  COUPLING  COEFFICIENTS 
FOR  THE  REFERENCE  MTERIAL  TEST  CONFIGURATIONS 


Energy  Coupling 

Coefficients  Thickness 

Configuration  (cal/g) / (cal/ cm^)  (cm) 


Sample 

Reference 

Sample 

^s 

Reference 

^r 

Sample 

Reference 

A1 

Qtz 

0.31 

0.33 

0.25 

0.26 

Ge 

Qtz 

0.32 

0.31 

0.20 

0.26 

Si 

Qtz 

0.34 

0.33 

0.30 

0.26 

Qtz 

Si 

0.34 

0.33 

0.26 

0.30 

Ge 

Si 

0.32 

0.31 

0.20 

0.30 
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TABLE  3 

REFERENCE  MATERIAL  PARAMETERS  USED  IN  THE  DATA  ANALYSIS 


Material 

Density 

(g/cm^) 

Wave 

Speed 

(cm/ys) 

Altminum  (6061-T6) 

2.70 

0.64 

Quartz  (X-cut) 

2.65 

0.572 

Germanium  (111) 

5.32 

0.550 

Silicon  (polycrystal) 

2.33 

0.593 

TABLE  4.  REFEREKCE  MATERIAL  DATA  S 


287  Ge/Si  0.0300  103  5.0  0.50  0.0560  106  19.1  0.81 


least-squares  straight  line  fits  to  the  pressure-energy  data  shown  in 
Figure  21;  they  are  normalized  to  be  consistent  with  the  value  Y = 2.1 
for  aluminm,  which  is  used  in  determining  the  deposited  energies.  The 
resulting  Gruneisen  parameters  are  given  in  Table  5 and  are  compared 
with  calculated  (thermodynamic)  values  and  with  results  previously 
measured  at  low  dose  (<  10  cal/g). 

An  alternate  approach  with  which  the  Gruneisen  parameters  were 
computed  from  incident  fluence  determined  "statistically"  from  the  fluence 
measurements  which  were  interspersed  with  data  shots  provided  less  con- 
sistent results.  This  result  is  attributed  to:  (1)  lack  of  machine 
repeatability,  and  (2)  inaccuracies  in  the  fluence  measurement  process. 
Extrapolating  the  peripheral  fluence  calorimeter  data  to  obtain  fluence 
at  the  central  sample  location  also  provided  inconsistent  results.  These 
results  are  exemplified  by  the  fluence  data  for  a series  of  shots  presented 
in  Figure  22.  This  figure  presents  the  axial  fluence  as  determined  by 
(1)  direct  measurement  by  carbon  fluence  calorimeter,  (2)  by  extrapolation 
of  peripheral  fluence  calorimeters,  and  (3)  by  computation  from  the 
reference  material  response  and  the  Gruneisen  parameters  of  this  work 
given  in  Table  5.  All  measured  (absorbed)  fluence  data  are  corrected 
to  incident  fluence  as  described  in  Section  4.  These  fluence  results 
emphasize  the  importance  of  obtaining  a direct  fluence  measurement  via  the 
reference  material  response,  on  each  specimen  test. 
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TABLE  5 


MEASURED  AMD  REFERENCED  GRUNEISEN  PARAMTERS 


Material 

'i  From 

This 

Work 

Y Calculated 
from  . 

Thermodynamic  Data” 

Y From 
Other  E-Beam 
Measurements' 

Aluminum 

2.10 

2.17" 

2.11^ 

2.14  - 2.19 

2.078 

Quartz 

0.65 

0.74® 

0.67^ 

0.65^ 

Germanium 

0.76 

0.72^ 

0.75® 

till] 

0.80" 

Silicon* 

0.50 

0.57^^ 

0,46j 

(polycrystal) 

0.43® 

0.42*^ 

a 

b 

c 

d 

e 

f 

g 

h 

1 

J 

k 

1 

m 

n 

ft 


Relative  to  Y.,  = 2.10 
A1 

Y = e/pKCp 

All  measurements  at  relatively  low  dose»  < 10  cal/g 

Reference  11 

Reference  3 

Reference  12 

Reference  13 

Relative  to  « 2.10»  Reference  3 
Reference  lA 

Relative  to  Y^  ■ 2.10»  Reference  3,  [111]  Si* 

Reference  13,  [111]  Si* 

Reference  IS 

[111]  Si*,  Reference  15 

Reference  16 

For  crystals  of  cubic  symmetry  such  as  SI  and  Ge,  the  same  Y 
is  obtained  for  a single  crystal  and  a polycrystal. 


47 


CHAPTER  7 
CONCLUSIONS 

1.  The  pressure-energy  relations  have  been  measured  to  100  cal/g 
for  four  reference  materials  (Al,  Qtz,  Si,  Ge).  Data  at  greater  doses 
would  be  of  value  in  extending  the  range  of  usefulness  of  these  reference 
materials. 

2.  The  four  reference  materials  can  be  treated  as  linear  and 
elastic  in  this  range  of  thermomechanical  response.  Measured  Gruneisen 
values  agreed  to  i 10  percent  with  calculated  values  using  thermodynamic 
properties.  The  following  Gruneisen  parameters  were  established: 

Aluminum  (6061-T6),  Y = 2.10 
Germanium  [ill],  Y » 0.65 
Quartz  (X-cut),  Y = 0.76 
Silicon  (polycrystal),  Y = 0.50 

3.  The  concept  of  using  the  reference  material  to  provide  an 
In-situ  measurement  of  dose  (fluence)  was  proven,  and  it  was  shown  to  be 
necessary  because  of  machine  non-repeatability. 

4.  The  concept  of  using  pressure  generating  reference  materials 
to  provide  more  accurate  measurement  of  initial  sample  stress  without 
heavy  dependence  on  the  sample  adlabat  was  demonstrated.  Part  II  of  this 
report  describes  the  application  of  this  concept  to  determine  ttie  Grimeisen 
parameters  of  porous  metals  and  cermets. 
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